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Spin ensemble based hybrid quantum systems suffer from a significant degree of decoherence 
resulting from the inhomogeneous broadening of the spin transition frequencies in the ensemble. We 
demonstrate that this strongly restrictive drawback can be overcome simply by burning two narrow 
spectral holes in the spin spectral density at judiciously chosen frequencies. Using this procedure 
we find an increase of the coherence time by more than an order of magnitude as compared to the 
case without hole burning. Our findings pave the way for the practical use of these hybrid quantum 
systems for the processing of quantum information. 


PACS numbers: 42.50.Pq, 42.50.Ct, 42.50.Gy, 32.30.-r 

Hybrid quantum circuits that conflate the advantages 
of different physical systems to achieve new device func¬ 
tionalities have recently shifted to the center of atten¬ 
tion [1]. This is largely because a new generation of 
experiments [MU lends encouraging plausibility to the 
vision of using such hybrid device concepts to reliably 
store and manipulate quantum information |12H17j . In 
particular, the recent achievements in strongly coupling 
large spin ensembles to superconducting microwave cav¬ 
ities [3® HH hold promise for combining many of the 
advantageous features of microwave technology with the 
long spin coherence times found, e.g., in crystallographic 
defects of diamond. 

Whereas the collective coupling to a whole ensemble 
of spins is the key to reach the strong-coupling limit, the 
ensemble generally comes with the downside of being in- 
homogeneously broadened, i.e., the transition frequencies 
between different spin levels are slightly different for each 



FIG. 1: (Color online) Sketch of the studied hybrid quantum 
system: a synthetic diamond (black) containing a spin en¬ 
semble (red arrows) coupled to a transmission-line resonator 
(curved gray line) confining the electromagnetic field to a 
small volume. 
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spin. As it turns out, the decoherence resultant from this 
broadening is currently the major bottleneck for the pro¬ 
cessing of quantum information in these hybrid quantum 
systems. First attempts at resolving this problem have 
meanwhile been put forward: On the one hand, it was 
shown that the decoherence is naturally suppressed for 
very strong coupling when the spectral spin distribution 
realized by the ensemble falls off sufficiently fast in its 
tails. Signatures of this so-called “cavity protection ef¬ 
fect” [nms], have meanwhile also been observed exper¬ 
imentally [niEo]. To fully bring to bear the potential 
of this effect requires, however, to go to very high val¬ 
ues of the coupling strength, which are presently difficult 
to reach experimentally. On the other hand, sophisti¬ 
cated concepts for the spectral engineering of the spin 
density profile have been proposed Eiiiia. These con¬ 
cepts rely, however, on a strong modification of the intrin¬ 
sically predefined density profile that is again very chal¬ 
lenging to implement experimentally. In this Letter, we 
present a method that circumvents the problems of both 
approaches by building on a very elementary concept that 
requires only a considerably reduced experimental effort. 
Specifically, we demonstrate that the burning of two ju¬ 
diciously placed spectral holes in the spin distribution 
suffices to drastically increase the coherence properties 
of the hybrid spin-cavity system. From the viewpoint of 
quantum control our approach constitutes a new and effi¬ 
cient strategy to stabilize Rabi oscillations in the strong¬ 
coupling limit of cavity QED [53H1S]- Suppressing the 
detrimental influence of inhomogeneous broadening, as 
suggested in our work, could also prove to be a key ele¬ 
ment for the realization of ultra-narrow linewidth lasers 

[13 HZ]- 

To connect our theoretical work directly with the ex¬ 
periment we will study in the following the recently 
implemented case of a superconducting microwave res¬ 
onator strongly coupled to an ensemble of negatively 
charged nitrogen-vacancy centers in a diamond (see 
Fig. 0 [3 E im [20] . Our starting point is the Tavis- 
Cummings Hamiltonian {h = 1) [28], which describes the 
dynamics of a single-mode cavity coupled to a spin en- 
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semble in the dipole and rotating-wave approximation, 

^ N N 

n = ujca'^a + 2 X! ^ a'*' - a] - 

j 3 

i [77(t)a^e-*“* - 7?(t)*ae‘‘"‘] . (1) 

Here cr^, cr“, cr| are the Pauli operators associated with 
the individual spins of frequency u)j. Each spin is coupled 
with a strength gj to the single cavity mode of frequency 
Wc, in which photons are created and annihilated through 
the operators and a. The probing electromagnetic 
field injected into the cavity is characterized by its carrier 
frequency ui and by the amplitude ri{t). 

Next, we derive the semiclassical equations of motion 
using the Holstein-Primakoff-approximation |29j (imply¬ 
ing that the condition (cr|) « —1 always holds), the 
rotating-wave approximation and neglecting the dipole- 
dipole interaction between spins. With these simplifica¬ 
tions, which are well justified for the experiments HHEo] 
operating at low input powers of an incoming signal, the 
equations for A{t) = {a{t)) and Bj{t) = {cr~ {t)) acquire 
the following form (in the w-rotating frame), 

A{t) = - [k -h i{ujc - w)] A{t) + - r]{t), (2a) 
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Bj{t) = - [7 + - w)] Bj{t) - gjA{t), (2b) 

where k, 7 are the dissipative cavity and spin losses, re¬ 
spectively. 

Large spin ensembles (TV ~ 10^^ in [TTl [^) are 
best described by the continuum limit of the normal¬ 
ized spectral density p(w) = 9~ Here 

H is effective coupling strength which 

is enhanced by a factor of \/N as compared to a single 
coupling strength, gj , so that H can reach the values nec¬ 
essary for the realization of the strong coupling regime. 
The inhomogeneous broadening of the spin frequencies 
uij and coupling strengths gj then lead to a finite-width 
distribution p{uj) centered around a certain mean fre¬ 
quency (Ms- The specific shape of this spectral density 
p{uj) can typically be determined by a careful compari¬ 
son with the experiment based on stationary [ 3 ] or dy¬ 
namical m transmission measurements. In the follow¬ 
ing we will use the same parameters as in [niEoj taking 
a g-Gaussian distribution [ 3 ] for p{uj) centered around 
uJsI‘2.'k = 2.6915 GHz, a full-width at half-maximum of 
= 9.44 MHz and a g-parameter of 1.39. The 
cavity decay rate, K/27r = 0.4 MHz (half-width at half¬ 
maximum) and the coupling strength H/27r = 8.56 MHz. 

The starting point for our strategy is the insight that 
the non-Markovian dynamics of the spin system, which 
is described by p(a;) and strongly coupled to the cavity 
mode, can be accurately modeled by an integral Volterra 
equation for the cavity amplitude A{t) (see Eq. ([ 5 ]) be¬ 
low and mi US]). The latter includes a memory-kernel, 
which is responsible for the non-Markovian feedback of 


the spin ensemble on the cavity, so that the cavity ampli¬ 
tude at time t depends on all previous events t < t. By 
performing a Laplace transform of this Volterra equation 
[ 20 ] or by carrying out a stationary transmission anal¬ 
ysis mi [IS] , the total rate of decoherence turns out to 
be r « /c -|- Tril'^p{uis ± H) in the limit of large coupling 
strengths, H > T and 7 —)■ 0. The value of T is thus 
determined by the spin density p(a;), evaluated close to 
the maxima of the two polaritonic peaks, a; = Wg ± H, 
split by the Rabi frequency flu « 2fl due to strong cou¬ 
pling. Our approach is now to take this relation literally, 
which is tantamount to saying that the decoherence rate 
r can be strongly suppressed by burning two spectral 
holes into the spin distribution p{u}) right at these two 
positions, ujh = tOs A: ft, such that p{iOh) = 0. The width 
of the holes A/i should be very small, such as to remove 
only a negligible fraction of the spins by the hole burn¬ 
ing. On the other hand, A/j is limited from below by the 
spin dissipation rate, A/j > 7- 

To demonstrate the efficiency of our approach ex¬ 
plicitly, we first perform a stationary analysis [H(t) = 
Bk{t) = 0] of the transmission T(w) through the mi¬ 
crowave resonator as a function of the probing frequency 
w. This quantity, which is directly accessible in the ex¬ 
periment mi ED], provides direct access to the occupation 
amplitude of the cavity [T{oj) oc H(a;)]. Assuming 7 —>■ 0, 
the transmission T(uj) acquires the following form. 


T{u;) 


in 

U! — UJc — fi^S{ui) + i[K + Trfl'^p{ui)] 


( 3 ) 


This expression is normalized such as to reach the maxi¬ 
mum possible value max(|T(w)|) = 1 for suitably chosen 
Lo, K, and p{uj). The real function S(uj) is the nonlinear 
Lamb shift ED] defined as 


S(u;) = V 



dujp(bj) 
lo—Cj ’ 


( 4 ) 


where V stands for the Cauchy principal value. In the 
reference case taken from the experiment [TTl [ 20 ] , p{ijj) 
has no holes, see Fig.j^a), and the transmission |T(w)p 
displays the well-resolved double-peak structure typical 
for the strong-coupling regime, see Fig. [^b). If we now 
burn two narrow holes into the spin density at the rele¬ 
vant positions = Wg ± fl, see Fig. [^d), and reevalu¬ 
ate |T(a;)p we observe a more than fiftyfold increase in 
the corresponding transmission peak values, see Fig.j^e). 
This dramatic change is all the more surprising consider¬ 
ing that the relative number of spins removed from p{u}) 
through the hole burning is less than 3%. 

To understand this behavior it is best to analyze the 
real and imaginary parts of the denominator of T{uj), 
see Eq. (§. For the observed transmission resonances 
at Cli = Wr with a maximum value of T{ujr) = 1 to 
occur, two conditions are satisfied simultaneously: (i) 
{ojr — uJc)/fl^ = < 5 ( 0 ;^) and (ii) p{u:r) = 0. Consider first 
condition (i): In the reference case without holes, see 
Fig-K c), the nonlinear Lamb shift i5(a;) displays rather 
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FIG. 2: (Color online) Comparison of the cavity coupled to 
the inhomogeneously broadened spin ensemble without and 
with hole burning in the spin density profile (left and right 
panels, respectively). Top row: The q-Gaussian spin density 
distribution, p(a)), without and with hole burning at luh = 
Both holes are of equal width, A(i/27r = 0.7 MHz, and 
have a Fermi-Dirac profile. Middle row: Transmission T{oj) 
without and with hole burning in p{u}) (note different y-axes 
scale). Bottom row: The corresponding nonlinear Lamb shift 
5(aj). Filled circles label resonance values ujr of the transmis¬ 
sion T{uj) occurring at the intersections between the Lamb 
shift 5(t<j) and the dashed line {uj — cUc)/!!^. At empty circles 
such intersections are non-resonant (see text). 


smooth variations in the vicinity of the resonant frequen¬ 
cies ujr, determined by the intersection of (5(0;) and a 
straight line (o; — ■ In contrast, for the case with 

hole burning, see Fig. if), (5(a;) exhibits rapid variations 
around the two resonance points within a very narrow 
spectral interval. As a consequence, the resultant trans¬ 
mission peaks become substantially sharper. Due to the 
second condition (ii) they also dramatically increase in 
height. Note, that no resonance occurs at o; = a;c be¬ 
cause p{uj) has a maximum at this point and condition 
(ii) is strongly violated, see Fig. [^c),(f). A close ex¬ 
amination of the structure of T(a;j shows, furthermore, 
that the narrow transmission peaks resultant from the 
hole burning do not replace the broad polaritonic peaks 
present in the reference case, but rather get to sit on top 
of them, see Fig.|^a). As will be seen below, the different 
resonance widths in T{uj) set two different time scales in 
the dynamics with, in particular, the sharp peaks in the 
transmission giving rise to an asymptotically slowly de¬ 
caying dynamics with a strongly suppressed decoherence. 

To explore whether the narrow holes we burnt into the 
spectral spin distribution at = WgiD have, indeed, the 
optimal location, we now also test all possible other hole 
positions symmetrically placed around the maximum of 
piuj) at u; = Ws. In Fig. [^b) we present the numeri¬ 
cal results for T{u}) as a function of the probe frequency 
w and of different hole locations iOh = ± d) scanned 



tj/2n [GHz] 



2.6795 2.6835 2.6875 2.6915 
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FIG. 3: (Color online) Transmission through the cavity T{uj) 
versus probe frequency uj for different locations of the holes, 
ujh, in the spin density profile, p{uj) (the width of the holes is 
A(i/27r = 0.7MHz). (a) Red (gray) curve: |T(aj)|^ in lin-log 
scale versus lo for ujh = uJs H. Black curve: Transmission 
in the absence of hole burning, (b) Yellow (light gray) areas 
mark the most pronounced peaks in |r(a;)|^ in the presence of 
hole burning. Blue (gray) areas stand for the secondary po¬ 
laritonic peaks which stem from the case without hole burn¬ 
ing. Dashed arrows designate the distance Q.r between po¬ 
laritonic peaks. The white vertical cut corresponds to the 
transmission shown in (a). 


between d; = 0 and d) = 16 MHz: While for large hole 
spacings (d) > 11.5 MHz) the effect of holes is negligi¬ 
ble, in the interval 0.8MHz<d,< 11.5 MHz we always 
find two sharp peaks superimposed on the two polari¬ 
tonic peaks approximately at the hole positions. Close to 
ujh = Ws ± these peaks are most pronounced and reach 
unity. In the limit when the holes are burnt very close to¬ 
gether (d) ;< 0.8 MHz) the sharp peaks merge into a single 
one, located directly at the central frequency iOs with a 
transmission maximum reaching again unity in the limit 
of (jJh —t Ws [see the yellow cusp in Fig. [^b)]. Using the 
symmetry of p{ui) with respect to cus, this behavior can 
also be proven analytically (not shown). To check the ro¬ 
bustness of our method we also tested different functional 
forms for the hole prohles (Fermi-Dirac, g-Gaussian and 
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FIG. 4: (Color online) (a),(b): Decay of the cavity occu¬ 
pation N{t) = from the initial state, for 

which a single photon with frequency oJc resides in the cav¬ 
ity and all spins are unexcited. The asymptotic decay Ge”’"* 
with and without hole burning (see red lines) is determined 
by the constants T/2tt = 3 MHz in (a) and a drastically re¬ 
duced r = 0.42k = 27r-0.17 MHz in (b). (c),(d): Dynamics of 
|H(t)|^ under the action of eleven successive rectangular mi¬ 
crowave pulses of duration corresponding to the Rabi period, 
T = 2n/QR = 52 ns, phase-switched by tt (every second pulse 
is shown as a vertical gray bar). Also here the asymptotic 
decay is much slower due to the presence of the holes. In all 
panels the holes in p(cu) have a width A^/27r = 1.4MHz and 
are burnt at t = 0 at wji = ojs ± D. 


rectangular distributions) and found qualitatively similar 
results to the Fermi-Dirac form employed for all of the 
above figures. 

To reach our ultimate goal of understanding the influ¬ 
ence of the spectral hole burning on the resultant dynam¬ 
ics, we now study the time evolution of A(t) explicitly 
for the resonant case u = uJc = <^s- The expression for 
the corresponding Volterra equation can be derived from 
(see [ 10 ] for details), 

A{t) = —KA{t) — (5) 

t 

j dujp{w) J - r;(t). 

0 

To prove that our predictions are valid not only in the 
semiclassical but also in the quantum case, we consider 
the case when all spins are initially in the ground state 
and the cavity mode a contains initially a single pho¬ 
ton, |I,4,). It can be shown that the probability for a 
photon to reside in the cavity at time t > 0, N{t) = 
(l,4,|a'^(t)a(t)|l,4,), reduces to N{t) = |(0J,|a(t)|l,4,)p = 
IA(t)p, where A{t) is the solution of Eq. with the ini¬ 
tial condition A{t = 0) = I (external drive ri{t) = 0). For 
the case without hole burning this solution is represented 
by the damped Rabi oscillations [see Fig. |^a)] found al¬ 
ready previously [ni ED]. By burning narrow holes in 
p{uj) at ujh = ujg ± n (immediately before t = 0 ), we ob¬ 
serve very similar transient dynamics, which is followed. 


Eqs. {2m 2b 


however, by a crossover to Rabi oscillations with a much 
slower asymptotic decay [see Fig. |^b)]. Quite remark¬ 
ably, the total decay rate F in this asymptotic time limit 
can even be substantially smaller than the cavity decay 
rate k alone. This is all the more surprizing since n was 
identified as the minimally reachable value for F in re¬ 
cent studies on the cavity protection effect [nnnEDi. 
Apparently a new type of physics is at work here: Al¬ 
though the system is in the strong coupling regime, the 
two spectral holes slow down the leakage of the energy 
stored in the spin ensemble back into the cavity. In par¬ 
ticular, when being even slower than the inverse of the 
cavity decay rate k, this sets a new global time scale 
for F, corresponding to the width of the sharp resonance 
peaks which we identified before in Fig. ia). From the 
mathematical point of view such a slow asymptotic be¬ 
havior is associated with the contribution of two poles 
in the Laplace transform of Eq. (H [ 20 ] . which appear 
when the holes in p{ijj) reach a critical depth. The pole 
contributions also stabilize the long-time behavior when 
the holes are shifted away from the polaritonic peaks [see 
Fig. I^b)], but the optimal hole positions remain close to 
the polaritonic peaks. Note that despite the consider¬ 
able photon loss {N{t) <C I) for long times the phase 
coherence is very well preserved here, a clear signature of 
which is the stable form of the Rabi oscillations. In this 
way a high “visibility” can be achieved, as required for 
the efficient processing of quantum information m- 

To demonstrate the efficiency of the hole burning ef¬ 
fect also for quantum control schemes, we pump the cav¬ 
ity by a sequence of tt phase-switched rectangular pulses, 
each with a duration corresponding to the Rabi period, 
T = 2'k/^1{i and a carrier frequency oj = cOc = ojg- As 
shown in m, this procedure is very well suited to feed 
energy into the strongly coupled cavity-spin system, lead¬ 
ing to giant oscillations of both spin and cavity amplitude 
[see left parts of Fig.|^c,d)]. Not only do we observe that 
these driven oscillations are more pronounced when burn¬ 
ing holes at Wft, = Wg ± n, but we find, in particular, that 
the Rabi relaxation oscillations setting in after switch¬ 
ing off the driving field are dramatically more long-lived 
than in the case without holes [compare right parts of 
Fig. ,d)]. These results confirm the robustness as well 
as the general applicability of our approach for various 
coherent-control schemes in the strong-coupling regime 
of cavity QED. 

In summary, we present an efficient method to sup¬ 
press the decoherence in a single-mode cavity strongly 
coupled to an inhomogeneously broadened spin ensem¬ 
ble. By burning narrow spectral holes in the spin den¬ 
sity at judiciously chosen positions the total decay rate 
is dramatically decreased to values that may even lie be¬ 
low the dissipation rate of the bare cavity. Experimen¬ 
tally, our approach can be implemented by exposing the 
cavity to high-intensity microwave signals with spectral 
components near the desired hole positions. Due to the 
strong driving the spins at these frequencies will equally 
populate their ground and excited state and will thus be 
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effectively removed from the coupling process with the 
cavity. 
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